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We study the correlation between the neutrino oscillation and the production/decay of the lightest 
supersymmetric particle(LSP) in the context of the minimal supergravity model(mSUGRA) without 
R-parity. We show how the neutrino masses and mixing which are consistent with the recent neutrino 
data can be obtained in this model, and describe how to probe the model by observing the LSP 
decay at collider experiments. It is shown that the generic 1-loop contributions to neutrino masses 
are too small to account for the solar and atmospheric neutrino oscillations and thus some fine 
cancellation in tree-level contribution is required to severely constrain the viable parameter space. 
In most parameter space of mSUGRA, the neutralino or the stau is LSP. Examining both cases, we 
find that there is a simple correlation between the neutrino mixing angles and the LSP branching 
ratios in the small tan /3 region so that the model can be clearly tested in the future colliders. In 
the large tan/3 region, such a correlation is obscured by the large tau Yukawa contribution which 
makes it nontrivial to test the model. 

PACS numbers: 12.60.-i, 14.60.St 



I. INTRODUCTION 



The recent progresses in neutrino experiments has led 
us to convince the existence of neutrino masses and fla- 
vor mixing. 0. In this regard, of primary interest 
are to look for New Physics candidates which are in- 
herited not only with a natural mechanism to generate 
the observed neutrino mass matrix and but also with 
some other predictions that can be tested in the near fu- 
ture. One of the best-motivated models endowing such 
a property would be the Minimal Supersymmetric Stan- 
dard Model (MSSM) with lepton number violation via 
R-parity breaking terms 01 as such models may 

produce lepton flavour violating signals predicted by the 
observed neutrino mixing and thus can be tested in fu- 
ture collider experiments . The particle spectrum in 
MSSM depends on the supersymmetry breaking mech- 
anism. One of the most popular scenarios for super- 
symmetry breaking is the minimal supergravity scenario 
(mSUGRA) [3 , which assumes a universal gaugino mass 
Mi/2, a universal scalar mass mo, a universal trilinear 
coupling Aq and Bo at Mqut scale. This framework 
can successfully evade the appearance of dangerous flavor 
changing neutral current(FCNC) and is highly predictive 
because there exist only 5 independent parameters in the 
model. In this work, we do complete phenomenological 
analysis on the masses and mixing of neutrinos and lep- 
ton number violating signature of LSP decays in the con- 
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text of mSUGRA with trilinear R-parity violation. The 
neutrino mass matrix is calculated up to the one-loop 
order, and it is studied how the predictions concerned 
with neutrino parameters in mSUGRA can be probed 
from the collider signals. In almost all parameter space 
of mSUGRA, the LSP is predicted to be a neutralino or 
a stau. In both cases, we calculate the production cross 
section, decay rate, and branching ratios of the LSP to 
investigate the correlations between neutrino oscillation 
parameters and collider signatures from the various chan- 
nels of LSP decay. 

II. NEUTRINO MASS MATRIX AND FLAVOR 
STRUCTURE OF TRILINEAR COUPLINGS IN 
MSUGRA 

Let us begin by writing the superpotential in the basis 
where the bilinear term LiH\2 is rotated away : 

W = ^H 1 H 2 + htL l H 1 E c l +hjQ l H 1 D c l +h l tQ l H 1 Ut (1) 



W = \UL Z E% + KL^Dl, (2) 

where Wo is R-parity conserving part, whereas W is R- 
parity violating part, and we assume the dominance of 
AJ = A^ 33 and \ = A^33 [i = 1, 2, 3) over the other trilin- 
ear couplings of lower generations and the universality at 
the GUT scale is imposed. In this framework, there are 
10 free parameters of the model, i.e. five conventional 
ones plus five R-parity violating ones; 

m , A), My 3 , tp, sign(» ; Ai 23 , Ai, 2 
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where tp = tan/3 is the ratio of two Higgs VEVs. As 
is well known, non-universality is developed at the weak 
scale via RGE evolution. The RGE in this basis can be 
found in Ref. Non- vanishing soft terms are given by 



V 



l LiH! LiH 1 



-BH 1 H 2 + h.c. 



(3) 



where Bi is the dimension- two soft parameter. We note 
that non-trivial VEVs for sneutrinos can be generated in 
this case. Including 1-loop contributions Vt oop , we obtain 
the following relation which comes from the minimization 
condition for the scalar potential, 
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where ^ = (i>i) / '{H®), and E^.' 2 '' are 1-loop contributions 



generated from Vi oop . 
presented in Ref. @,| 



The explicit forms of si 1 ' 2 *' are 
. Introducing another variables, 



m = 6 

we can obtain the relation, 
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whose non-zero values are due to the neutral scalar loops. 
Then, the tree-level mass is presented in terms of £j; 
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where F N = M 1 M 2 /{c 2 v M x + s^Ma) + M|c 2/3 //x. In- 
cluding 1-loop corrections, the neutrino mass matrix is 
written as 



-f^G cos 2 /3 
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where Hij denotes the 1-loop contribution of the neutrino 
self energy, and 
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Exact expressions of the 1-loop contributions II' s are 
presented in Ref. @- We note that the conventional 
tau-stau, bottom-sbottom loops and neutral scalar (sneu- 
trino/neutral Higgs boson) loops are essential to achieve 



FIG. 1: Plots tan/3 vs. mass square ratio for general points, 
the region between the two straight lines is the allowed one 
by neutrino data 



a realistic neutrino mass matrix for tan/3 < 30, and they 
are given by 



M loop = 
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It is found that the charged scalar (slepton/charged Higgs 
boson) loops are important only for tan/3 > 30 as ob- 
served in Ref. whereas the neutral scalar loops are 
important for all tan /3 in our case. We have observed 
that the tree level masses are much larger than the 1- 
loop contributions in mSUGRA, which makes it difficult 
to accommodate the solar and atmospheric neutrino os- 
cillations. Thus, we need some cancellation in so that 
the tree mass becomes comparable to the 1-loop contri- 
butions. Then, we can obtain a desirable neutrino mass 
matrix the by combining the tree and 1-loop masses ap- 
propriately, but it is possible in some fine-tuned param- 
eter space, as will be shown later. In this case, however, 
we lose nice predictability of atmospheric and solar neu- 
trino mixing angles measurable in colliders. 



III. NUMERICAL RESULTS: FITTING THE 
NEUTRINO DATA 

First of all, we obtain the parameter sets consistent 
with the solar and atmospheric neutrino data by scanning 
the input parameters within the ranges, 

lOOGeV < m < lOOOGeV, (12) 

lOOGeV < Mi /2 < lOOOGeV, (13) 

OGeV < A < 700GeV, (14) 

2 < tan/3 < 43. (15) 

The ranges for R— parity violating parameters we scan 
are given 



4 x 10~ 6 < |Ai| < 6 x 10" 
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FIG. 2: Atmospheric neutrino mixing angle vs. ^3/^3! and 



A2/A3 



for general points 
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FIG. 3: The solar mixing angle and Ai / A2 ratio for solution points, Left : small tan (3, ~ 3—15 Right : for large tan /3, ~ 30 — 40 
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We set the signs of Mi/ 2 and A arbitrary, but find 
that most of the allowed parameter space corresponds 
to the case that both signs are positive. We also find 
that there are strong correlations between atmospheric 
neutrino mixing angle and A/s. In order to account for 
the large atmospheric angle and small Chooz angle, we 
should take X 1 <C A 2 ~ A 3 . We have confirmed this by 
numerical calculations. In mSUGRA scenario, the tree 
level value of neutrino mass is much larger the 1-loop 
contributions in general, whereas the allowed parameter 
sets from the neutrino data are not tree-dominant at all. 
So we should make the tree mass terms comparable to 
the 1-loop contributions, and it can be possible in the 
case that some cancellation in xii occurs In doing so, 
however, we lose beautiful predictability of atmospheric 
and solar neutrino mixing angles measurable in colliders. 
FIG.l shows scattered plot of the ratio Am^/Am^ 
vs. tan/3 in the randomly generated parameter space. 
The points in the region between the two dotted lines 
give the right value of the mass square ratio. Only small 
set of scattered points lie in that region. We find that 



only small number of points among them are consistent 
with bi-large mixing. 

We note that the relative sizes of R-parity violating 
parameters are correlated with neutrino mixing pattern. 
But due to the cancellation in xii, it turns out that such 
a correlation is rather weak. In FIG. 2, we see some cor- 
relations between the atmospheric neutrino mixing an- 
gle and A 2 /A 3 for general parameter sets although 
they are not very strong because many data points do not 
lead to tree-dominant neutrino masses. Even though the 
correlations look so weak, we can obtain some constraints 
on |A 2 /A' 3 | from the neutrino data 0.4 < |A 2 /A 3 | < 2.5 
for small tan (3 and 0.3 < |A 2 /A 3 | < 3.3 for large tan/3 In 
FIG. 3, we plot the solar mixing angle vs. A1/A2 ratios 
for small tan/3 and large tan/3 cases, respectively. We 
see that there are strong correlations between the solar 
mixing angle and Ai/A 2 ratios. Similar to the above case, 
we can get some constraints 0.3 < IA2/A3I < 1.6 for small 
tan/3 and 0.2 < |A 2 /A 3 | < 5 for large tan/3. Collecting 
these constraints, the solutions are found for the ranges 
of 



Ai,2,A 2i3 = (0.1-2) x 10- 4 , 
A; < 2.5 x 10~ 5 . 



(21) 
(22) 



Let's summarize the results of this section. In general, 
in large portion of parameter space, the tree level con- 
tribution is much larger than those from 1-loop, so it is 
hard to get the right mass square ratio which is about 
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0(1O -2 ). The solutions are possible by suppressing the 
tree level contributions, which arises through the cancel- 
lation in £j. In this case, the correlations between £j and 
mixing angles become worse. In addition to the impor- 
tant 1 loop correction by stau and sbottom, the neutral 
scalar loop can give the important contribution to the 
neutrino masses, through the deviation of n from the 
direction of £j, which determine the tree level neutrino 
mass matrix. Instead of those, there are stronger corre- 
lations between the atmospheric mixing angle and A 2 /A 3 , 
the solar mixing angle and A1/A2. For small tan/3, it is 
possible to probe these in collider signal. It will be shown 
in the next section. 



IV. COLLIDER SIGNALS OF THE MODEL 

It is expected that the luminosity of the next linear col- 
lider can reach the 1000fb _1 /yr, and the center of mass 
energy over With such a capacity of the 

linear collider, it could be possible to prove the supersym- 
metric particles pair produced^] as weu as the decay of 
LSP inside the detector We present here the total 
cross section and decay rate for various parameter sets, 
and discuss the possibility to probe the structure of Im- 
parity violating parameters constrained by neutrio data. 
In the previous sections we got the parameter sets which 
are consistent with the known neutrino data. Those re- 
sults suffice to determine which is the LSP, neutralino 
or light stau. In TABLE I-III, we present the results of 
decay rate and branching ratios. The TABLE I shows 
the results for the stau LSP case with small tan (3. TA- 
BLE II and III correspond to the Neutralino LSP cases 
for small tan/3. The former is the case that only 3-body 
decay is permitted, whereas the latter is the case that 
both 2-body and 3-body decays are permitted. Since the 
decay lengths are smaller than a few cm for all cases, it 
is possible to detect their decay modes in the next linear 
collider. 

Firs of all, let's discuss the stau LSP cases. For small 
tan /3, fi ~ tr due to the small off-diagonal part. Then 
the light stau almost decays into leptons via XiLiL 3 E^ 
terms in the superpotential. In this case, the following 
relation holds, 

Br(ev) : Br(fj,v) : Br(rv) 

~ |Ax| 2 : |A 2 | 2 : |A!| 2 + |A 2 | 2 . (23) 

Thus, by observing the lepton branching ratios, one can 
measure the ratio of Ai and A 2 . If this value is out of the 
range constrained by neutrino data, given in the last sec- 
tion, we can exclude this model. If tan /3 become larger, 
the above relation does not hold any more because the 
Yukawa coupling plays an important role in this case. In 
other words, f\ cannot be considered as almost tr, but 
should be the combination of tl and tr and then the 
couplings are given by the mixture of Ai and h T . 

Next, let's consider the neutralino LSP cases. The 2- 



body decay rates are proportional to |£i| , 

T(yiZ) tx |£i| 2 , (24) 
T{hW) <x |6| 2 . (25) 

For 3-body decay modes, the lepton-quark-qurak branch- 
ing ratios are proportional to |£i| 2 since W boson ex- 
change diagrams are dominant. So we can have the fol- 
lowing relation, 

Br(ejj) : Br(fijj) : Br(rjj) = |£i| 2 : |£ 2 | 2 : I6| 2 - (26) 

Thus, we can obtain the information of £j ratio by mea- 
suring the branching ratios, but it is difficult to test 
mSUGRA scenario for massive neutrinos via collider sig- 
nals because the correlation between £j's and neutrino 
oscillation parameters is diminished as pointed out in the 
last section. Similar to the stau LSP case, for small tan /3, 
branching ratios provide Ai information on them; 

Br(ve ± T T ) : Br^^T*) : Br{vT^T*) 

~ |A!| 2 : |A 2 | 2 : |A!| 2 + |A 2 | 2 (27) 

Likewise, if we measure the above branching ratios, the 
models can be tested by comparing them with the range 
allowed by neutrino data, But since, for large tan/3, the 
terms generated from h T Lj,H\E 3 in the superpotential 
become large, and thus the above relation l|27|) breaks, 
it is impossible to test this scenario. The examples are 
presented in the TABLE II and III. In addition, we note 
that ?7i's are rather big in all cases, which implies that 
some cancellations in £j's occur, as expected before. 

V. CONCLUSION 

Based on the recent precise experimental results for 
neutrino oscillations, we have studied the minimal su- 
pergravity model with the lepton number violation via 
R-parity violating terms Q . Since it is impossible to ac- 
count for large mixing angle solution of solar neutrino 
problem in mSUGRA with only bilinear R-parity viola- 
tions, we should include the general trilinear R-parity 
violations. So there are additional 5 parameters Ai and 
Ai besides the 5 parameters which are from the ordinary 
minimal supergravity model. In future collider which can 
reach its CM ener gy o ver 1 TeV and luminosity up to 
1000 fb _1 in a vear|ll] [ l2 ^ . we expect that it is possible 
to probe mSUGRA without R-parity via the LSP de- 
cay. We have presented how neutrino mass matrix from 
mSUGRA with R-parity violation can be constructed at 
the 1-loop level. From our analysis, we found that most 
parameter space which is consistent with neutrino data 
leads to 1-loop contributions as large as comparable with 
tree level masses. We have found the strong correlation 
between solar neutrino mixing angle and Ai/A 2 , and also 
between atmospheric neutrino mixing angle and A 2 /A 3 . 
For searching the LSP decays, we have considered two 
cases in mSUGRA, neutralino LSP and scalar tau LSP, 
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setA : Stau LSP 




Ldll — <J.J_>J 




/; — —801 93 fieV 

fJi — OUl.^O VjrC V 










A' 


8.403 x 10" 6 


-7.663 x 10~ 5 


-6.792 x 10~ 6 




8.739 x 1(T 5 


-7.444 x 10~ 5 





& 


-1.003 x 10~ 6 


-3.642 x 10~ 6 


-4.401 x 10~ 6 




-9.897 x 10~ 7 


-3.045 x 10~ 6 


-3.728 x 10~ 6 


Vi 


1.011 x 10~ 6 


-1.313 x 10" 5 


-1.192 x 10" 5 


BR 


e 


A 4 


T 




28.93 % 


20.99 % 


50.02 % 


t b 




~ 0.066 % 






771^=278.59 GeV 


r = 


2.921 x 10~ 7 GeV 




^~e + e~ — >f i f * 


~ 1.450 xlO~ 2 (Pb), 


= 1 TeV 



(Am|i, Amii) = (2.50 x 10~ 3 , 1.13 x 10~ 4 ) eV 2 
(sin 2 2(9 atm , sin 2 2<9 so! , sin 2 29 chooz ) = (0.98, 0.77, 0.03) 



TABLE I: A trilinear model realizing the LMA solution for stau-LSP case. Here the couplings AJ and A, can be considered as 
input parameters defined at the weak scale. 





setB : Neutralino LSP 




tan/3 = 4.94 


sgn(^) = -1 


fj, = -200.46 GeV 




A = 38.93GeV 


m = 333.66 GeV 


Mi /2 = 160 GeV 


K 


-9.326 x 10" 9 


-7.811 x 10" 5 


-7.560 x 10 -5 


A, 


-5.628 x 10~ 5 


-7.345 x 10" 5 
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-1.234 x 10~ 6 


3.247 x 10~ 7 


1.880 x 10~ 6 




-1.211 x 10~ 6 


2.749 x 10~ 7 


1.831 x 10~ 6 


m 


-2.007 x 10" 6 


-1.169 x 10" 5 


-8.742 x 10" 6 


BR 


e 


A* 
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"33 




47.01 % 




Ifjj 


3.88 xlO' 2 % 


2.00 xl0~ 3 % 


8.87 xl0~ 2 % 




9.76 % 


16.60 % 


26.39 % 




771^0=59.37 GeV 


r = 


7.137 x 10" 15 GeV 






~ 4.897 xlO~ 2 (Pb), 


y/s = 1 TeV 



(Amii, Amli) = (2.51 x 10 -3 , 9.64 x 10 _s ) eV 2 
(sin 2 2<9 a tm, sin 2 2<9 soi , sin 2 26 chooz ) = (0.98, 0.99, 0.008) 



TABLE II: Neutralino-LSP case, only 3-body decays are possible 



and searched each parameter space regions by varying 
to an d Mi/2- For our purpose, we have calculated the 
production cross section for each cases. For stau LSP 
case, 2-body decay modes are dominant. Since stau LSP 
is almost tr for most parameter space, its branching ra- 
tios of the decay into top and bottom quark is negli- 
gibly small. Since stau LSP decay into liv takes place 
via Ai couplings, it is possible to obtain the informa- 
tion on Ai's from branching fractions. Wc found that 
Br{eu) : Br^v) : Br{ru) = |Ar| 2 : |A 2 | 2 : \\ x \ 2 + |A 2 | 2 
for not too large tan/3. For neutralino LSP cases, sim- 
ilar analysis has been done. Unlike the stau LSP cases, 



neutralino LSP can decay into 3-body final states as well 
as 2-body decay. ^From the branching fractions BrQijj) 
for various decay channels of neutralino LSP, we have 
determined the ratios of £j's, which in turn give the de- 
sirable tree level neutrino mass matrix. But, due to 
rather large loop contributions to neutrino masses the 
predictability for neutrino parameters is diminished. For 
small tan (3, we have obtained the Ai 's information which 
are strongly correlated with solar neutrino mixing angle 
from the branching fractions Briyl^ r^). But unluckily, 
since Yukawa couplings become larger for large tan/3, it 
is impossible to get information on Ai's. Consequently, 
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setC : Neutralino LSP 
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A, 


-7.966 x 10" 5 


-7.995 x 10 -5 





e? 


1.274 x (T 6 


2.896 x 10~ 6 


2.138 x 10~ 6 


& 


1.212 x 10" 6 


2.976 x 10~ 6 


2.409 x 10 -6 


T< 


4.038 x 1CT 7 


-1.305 x 10" 5 


-1.768 x 10" 5 


BR 


e 




T 


km 




14.77% 






6.36xl0" 2 % 


3.84 xl0 _1 % 


2.51 xl0~ x % 




20.96 % 


21.186 % 


42.09 % 




m.o =195.00 GeV 


r = 


4.744 x 10 -11 GeV 






~ 0.197 (Pb), 


= 1 TeV 



(Am§ x , At4) = (2.50 x 10~ 3 , 6.81 x 10~ 5 ) eV 2 
(sin 2 2(9 at m, sin 2 26» soi , sin 2 26 chooz ) = (0.95, 0.96, 0.007) 



TABLE III: Neutralino-LSP case, both 2- and 3-body decays are possible. 



for both stau LSP and neutralino LSP, it is possible to 
probe the relations between neutrino oscillation and Im- 
parity violating couplings from the branching fractions. 
But it becomes hard for large tan (3. Thus, any observa- 
tion of the sever deviation from these branching fractions 
in the future collider may exclude scenarios for neutrino 
masses in the framework of mSUGRA with appropriate 
R-parity violation. 
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